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Abstract 

In this paper we show that the quantum theory of chaos, based on 
the statistical theory of energy spectra, presents inconsistencies difficult 
to overcome. In classical mechanics a system described by an hamilto- 
nian H = Hi + H2 (decomposable) cannot be ergodic, because there are 
always two dependent integrals of motion besides the constant of energy. 
In quantum mechanics we prove the existence of decomposable systems 
H q = HI + H2 whose spacing distribution agrees with the Wigner law and 
we show that in general the spacing distribution of H q is not the Poisson 
law, even if it has often the same qualitative behaviour. We have found 
that the spacings of H q are among the solutions of a well defined class 
of homogeneous linear systems. We have obtained an explicit formula for 
the bases of the kernels of these systems, and a chain of inequalities which 
the coefficients of a generic linear combination of the basis vectors must 
satisfy so that the elements of a particular solution will be all positive, 
i.e. can be considered a set of spacings. 

1 Introduction 

When one talks about an ergodic mechanical system in classical mechanics, 
everybody refers to the same definitions and the same theorems. In quantum 
mechanics, unluckily, a universal agreement on quantum chaos doesn't exist. Up 
to now several criteria have been proposed for discriminating chaotic hamiltoni- 
ans from ordered ones, which either present limits and inconsistencies difficult 
to overcome, or need further evidence to have the same consistence of the clas- 
sical theory |g [k], 0, H, || |l4| ||, ||| . There is still a lot of confusion on the 



subject, so that one could even doubt the very existence of this vague concept. 
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The mathematical properties characterising ergodicity in quantum mechan- 
ics, if they exist at all, are still way ahead of us. Whatever they are, however, 
they ought to be intrinsic, i.e. independent of the corresponding classical sys- 
tems, which, in general, could even not exist. Notwithstanding this, it is rea- 
sonable that when there is a classical limit, an ergodicity criterion for quantum 



systems give answers which agree with the classical behaviour 12 . Since in clas- 



sical mechanics two non interacting hamiltonians cannot form a chaotic system, 
one should ask that this be true also for the quantum counterpart, whatever 
criterion he is using; otherwise serious doubts would rise about its validity. 

In the light of this simple test we have analysed the most popular criterion 
for quantum chaos, which is based on the study of the statistical properties 
of energy spectra j|, ||, [|: it is commonly believed that, if the distribution of 
spacings between neighbouring levels obeys Poisson law, the system is ordered; 
otherwise, if the they are distributed according to Wigner law, the system is 
chaotic. 

Clearly the analysis we want to perform is esentially algebraic. Let the 
eigenvalues of H q be 

E ij =El+E* '../ 1.2 //. (1) 

First we have studied the conditions which Eij must satisfy so that the 
algebraic system (Q) can be solved in the unknowns E\ and then, carrying 
out lengthy algebraic calculations, we have found that the spacings of H q are 
among the solutions of a well defined class of homogeneous linear systems. We 
have obtained an explicit formula for the bases of the kernels of these systems, 
and a chain of inequalities which the coefficients of a generic linear combination 
of the basis vectors must satisfy so that the elements of a particular solution 
will be all positive, i.e. can be considered a set of spacings. 

Solving these inequalities seems to be a formidable task, and we have not yet 
been able to find a way to do it. The solution of this problem is the starting point 
for classifying the weights of the possible histograms of spacings of decomposable 
systems. Numerically we have found histograms quite different from Poisson law 
and also some consistent with the Wigner distribution. We think, even if we 
cannot yet prove it, that their weight is not negligible. 



2 Conditions on the energy levels 

Eq. (|l|) in the previous section is a linear system of n 2 equations in the 2n 
unknowns E\ and E 2 , whose coefficient matrix has the following structure: 
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/ is the identity matrix, while S , ' m ) is a matrix whose elements are all except 
those in the ra — th column whose entries are all 1, 
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It is easily seen that the rank of the matrix 
dependent equations is n 2 — 2n + 1 = (n — l) 2 . 

It is important to examine carefully the properties of systems whose energy 



1Z\ , whose 



levels belong to the range of the linear operator A4 
matrix realization in the canonical basis is given by (^|). 

Let M Q be the n 2 by 2n + 1 matrix obtained adding to M a column vector 
E containing the energy levels of the system H q = Hi + H\. Reducing by rows 
the matrix (0), one obtains 
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is a n by n matrix with only one element different from zero in the first 

row: 

rr( m ) _ x 

Uij — Ol l+j-m- 

Repeating the same sequence of operations leading to matrix (||) on the matrix 
M Ql and imposing that it have the same rank of M, it can be shown that 
the elements of the vector E must satisfy the homogenous linear system whose 
coefficient matrix is the following: 
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where B is the matrix 
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The matrix C has (n — l) 2 rows and n 2 columns; it can be decomposed in a 
natural way in blocks of matrices of n — 1 rows and n columns. It is reduced by 
rows; therefore, its kernel has dimension 2n — 1 and coincides, by construction, 
with the range of M.. A basis of the kernel obtains by considering 2n — 1 
arbitrary columns of the matrix (||). The analitical expression of (|J) is 

Bki = $n — 5 k+ u, 

Cx^m — (Sifj, — 6a+i^i) Bu, 
A, k — 1, 2, . . . , n — 1, 
H,l = 1,2, ...,n. 

The indices A, \i denote the position of a block in the matrix C, while fc, I locate 
the single element. They are related to the usual indices of rows and columns 
r, s by the relations 

r = (A - l)(n - 1) + k, 
s = (/! — l)n + I. 



3 Connection between energy levels and spac- 
ings 

In the previous section we have treated the decomposability of a spectrum; yet 
the statistical analysis concerns the spacings. Our goal is to know when, given 
the number of elements in each cell of a histogram, we can find n — 1 spacings 
- compatible with the histogram — which may belong to a decomposable 
system. 

A set of n 2 — 1 spacings can be consistent with (n 2 — 1)!, generally distinct, 
spectra (just this particular feature might bring about doubts on the effective- 
ness of the spacing analysis as a tool for probing quantum systems). We can 
create a different system of levels, apart from an unimportant costant, for every 
single sequence of spacings; without loss of generality, we will set the arbi- 
trary constant, which coincides with the eigenvalue of the fundamental state, 
equal to zero in every physical system: in Hf, iJ| and in the total hamiltonian 
m = Hi + Hi 

With the above convention, a possible system of n 2 levels is obtained apply- 
ing the following n 2 by n 2 — 1 matrix to a vector of spacings: 
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Figure 1: Distribution of spacings of a decomposable sandwich system. 

For E = TA, where A is the spacing vector, the decomposability condition 
becomes 

CTA = 0. (6) 

T is not the most general matrix which may connect a set of spacings with 
the levels of a decomposable system; it follows from what has been stated above 
that, if one arbitrarily permutes the columns of (J5|), one obtains a new matrix 
T c which connects the same set of spacings with another system of levels among 
the possible (n 2 — 1)!. Furthermore by the structure of the matrix (Eh it is plain 
that the elements Ej obtained through T — or T c — are increasing with the 
index j. Given the structure of the matrix (g), and our convention, the elements 
of vector E are increasing with the index j if and only if the spectrum of the 
second hamiltonian is entirely contained within every spacing of the first one, 
i.e. if they look like a ladder "sandwich" . 

This is, of course, a very particular case, but it is indicative, we believe, of the 
possibilities that may occur: as a matter of fact we can always choose Hf and 
-fff so that the distribution of spacings of H q will be consistent with an arbitrary 
one: Poisson, Wigner or whatever we like. The histogram of spacings of H q is 
obtained, but for n elements, multiplying by n — 1 the number of spacings in 
the cells of the histogram of fff, an d it will have the same qualitative behaviour 
(see fig. 0). 

A matrix connecting a decomposable vector E, whose elements are not in- 
creasing, is obtained permuting the rows of (g), but, as we shall see, only the 
pemutations of a particular class are allowed for rows. 
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If, without loss of generality, we arrange the levels of Hf and iff in increasing 
order, the eigenvalues of the total system must obey the following inequalities: 

Eij+i > Eij, Ei + ij > Eij] (7) 

the double indices have the same meaning as in eq. (]]]). 

Let us take a decomposable system of levels, not necessarily sandwich, and 
arrange the eigenvalues in increasing order, 

E t > Ej if i > j; (8) 

such a vector belongs to the range K of Ai, or to the kernel of the matrix ([IJ), 
only if the system is sandwich: in all other cases an appropriate permutation 
exists such that the vector E a u\ G K. Inequality (|J) implies that the elements 
of the new vector must preserve the original order, i.e. 

E a{i) > E a[j) if i > j. (9) 

A level of the total system can be denoted in two ways: with a double index 
— Eij — which shows that it is the sum of the i — th level of the first system and 
the j — th level of the second one, or with the usual single index that locates 
the element in vector E. Therefore it is natural to single out the following 
association: 

(1.1) - 1 

(1.2) - 2 



(l,n) — > n 
~ n + l 
-> n + 2 

(10) 



(2,1) - 


n + l 


(2,2) -> 


n + 2 


(2,n) -> 


2n 



(n, 1) -> n(n - 1) + 1 
(n,2) -> n(n-l) + 2 



(n, n) — > n 2 

This map will work for sandwich systems only; for a generic decomposable 
system we have to replace the second column with an appropriate permutation. 
The homogeneous system (g) now becomes 

CT p A = 0, (11) 

where T p is obtained from matrix (^) by means of a suitable permutation of rows. 
In general with an arbitrary permutation the system (11) may not admit phys- 



ical solutions, i.e. the spacings could be negative, or, equivalently, the system 
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CE = doesn't have solutions verifying inequalities (^). The problem of locat- 
ing the allowed permutations is equivalent to finding a particular class of maps 

A [ J : N e i — > TV. 

It is useful to denote the range of A V by means of a n by n matrix A^ p \ 
From (|J) and (|l^) it follows that a permutation is allowed if and only if 

4O) ^ 4O) j a (p) ^ 4 GO. 

this subclass of allowed permutations in S^i 2 is made up by all and only those 
maps which can be allocated in a matrix in such a way that the elements of the 
rows are strictly increasing from left to right and the elements of the columns 
are strictly increasing from top to bottom. It is well known that the number of 
permutations having this property is the dimension of the irreducible represen- 
tation characterized by the corresponding square Young tableau. 

It follows from the above considerations that the decomposability of a non- 
negative vector A is expressed in the most general way by the homogeneous 
systems (pi]), in which matrix T p is obtained from T by means of a generic 
permutation of the columns and an allowed one of the rows. 

It is important to point out, for what will follow, that the allowed permu- 
tations leave always the first and the last element fixed, so the elements of the 
first row of T p are all zero. 

A detailed study of the kernels of all matrices CT p is necessary in order to 
classify all possible histograms of decomposable systems. 



4 The chain of inequalities 

The solution of the system (pT|), to make sense, must be made up of nonnegative 
elements. It follows that in our analysis we have to consider only the intersection 
of the kernel K p of the matrix D p — CT p with the positive hyperoctant of 
TZ^~°°. An arbitrary basis of a subspace is not, in general, formed by vectors 
whose elements are all nonnegative; moreover it is difficult to know a priori if 
a nonnegative vector exists in such a subspace, and, if it does, it is much more 
difficult to find it. Nevertheless it is easy to prove that if a vector of K p belongs 
to the positive hyperoctant (origin and coordinate hyperplanes excluded) then 
a whole basis made up of nonnegative vectors exists. 

Let us find a basis of K p . Let T be a triangular matrix obtained removing the 
first row from matrix (^), and let T~ 1 be the matrix obtained permuting the rows 
and the columns of T -1 by means of the permutations applied to the columns 
and to the rows of (||), respectively, to obtain T p . The allowed permutations 
contain n 2 elements, while the square matrices T^ 1 have dimension n 2 — 1. 
This however is not a problem because the first row has been removed from 
the matrix T, and the allowed permutations do not change the first element; 
therefore to obtain the right permutation to apply to the columns of T _1 it is 
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enough to remove the first element of the corresponding allowed permutation 
and to subtract one from the other elements. So the matrix T p coincides with 
T p where the first row has been removed. 

We have already noted that the kernel of (0) coincides with the range of (||) . 
A basis of this subspace is made up of 2n — 1 arbitrary columns of the matrix 
M. Let M be the matrix obtained from M removing the first row: the matrix 



Tpf^M 



(12) 



concides with M , except for the first row which now has every element equal 
to zero. Let N be the submatrix of M obtained removing the first and the 
(n + 1) — th column; the matrix 



W = f^N 



(13) 



is made up of columns which, by construction, are vectors belonging to K p ; 
since T~ x is obviously invertible, such vectors are linearly independent. We 
have excluded the first and the (n + 1) — th column of the matrix M, because 
their images through T p T~ x are the only columns of ( |l2| ) which differ from the 
corresponding columns of M, and therefore do not belong to the kernel of (|J). 
It is an easy exercise to prove that the vectors of K p so obtained are a maximal 
set of linearly independent vectors, i.e. that they are a basis. 

The kernels of matrices D p are subspaces changing for different choices of 
the permutations applied to the rows and to the columns of (||): what we are 
interested in is a non-empty intersection of K p with the positive hyperoctant, 
it should be clear by now, that a permutation of the rows of (||) can shift 
this intersection outside the positive hyperoctant, whereas a permutation of the 
columns can only shift it inside. 

Now we will find the conditions which must be satisfied by the coefficients 
of an arbitrary linear combination of the columns of ( |l3| ) so that such a vector 
will be nonnegative; obviously we will restrict our analysis to matrixes T p whose 
rows are changed by allowed permutations. 

The matrix T -1 has the following structure: 
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its analitical expression is the following 
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the analitical expressions of the columns of the matrix TV are the following 

n 

m i = ^2 S in(k-l)+r r = 1,2, . . .,n- 1, (15) 

n-1 

m n-i+r = ^2 § . nr+k i^l, 2 ,...,n 2 -l. (16) 

Let Tp 1 be obtained permuting the columns of T~ l by means of an allowed 
permutation; multiplying on the right by N we have 

f~ x N = f- 1 ^, (17) 

where N p i is obtained permuting the rows of N with the inverse permutation 
of that applied to the columns of T _1 . Now formulae ( [l5|) and (|l6| ) become 

n 

m i = E^" I (')»H)-h' r = l,2,...,n-l, (18) 
fe=i 

n-l 

m?- 1+r = X) <y «- 1 (0nr+* i = l J 2,...,n 2 -l; (19) 

a~ l {i) denotes the inverse permutation of the allowed one. Finally we obtain 
Wi r = 8 a -i^) nq+r — 5 a -\^x) np j rr i = l,2, ...,n — 1, (20) 

Wir+n-l = ^a-^i) nr+p ~ ^-'(i-l) nr+9 f = 1, 2, . . . , 71 — 1; (21) 

p, q, p and q are indices running from to n — 1. These expressions show that the 
second member of equation ( |20|) is different from zero if the index r coincides 
with the remainder of the division by n of a _1 (j) or a _1 (j — 1), while the 
second member of equation (^l| ) is different from zero if the index r is equal to 
the quotient of the same division. If we require that the individual elements of 
an arbitrary vector, linear combination of the columns of W , be strictly positive, 
we obtain the following chain of n 2 — 1 inequalities for the 2n — 2 coefficients of 
the linear combination: 

+ C Q a -Hi) > CR a -Hi-x) + C Q a -Hi-i) > ( 22 ) 

Ra- % {i) denotes the remainder of the division by n of a~ l {j), while 

its quotient. C s and C s are the coefficients of the columns of W, with the 

convention that they are zero when s is equal to zero. 



1 Because these are basis vectors, we find it convenient to use a notation in which the upper 
index denotes the k — th vector, and the lower index the i — th element of such a vector. 
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We can conclude that if the chain (|22|) admits solutions the kernel of the ma- 
trix Dp intersects the positive hyperoctant and the solutions of inequalities ( p2| ) 
characterise completely the intersection. 

The problem of classifying histograms compatible with such an intersection 
is still unsolved. In a histogram the choice of the width of the cells depends 
on the average number of elements per cell; in unfavourable cases, when the 
number of samples is small, few cells or few elements per cell yield a mean- 
ingless histogram, because a different choice of the cells changes completely its 
behaviour; when, on the contrary, the number of samples is large a variation of 
the width of the cells, in a reasonable interval, leaves its appearance unchanged. 
In what will follow we shall consider two histograms distinct if their cells have a 
different number of elements; clearly two different histograms in our convention 
may represent the same distribution, but, given the freedom in the choice of the 
cells and the algebraic nature of the problem, it is difficult to find a different ap- 
proach. Nevertheless if one could carry out this classification one would succeed 
in identifying all the distributions compatible with decomposable systems. 

In our analysis we will choose the unit of the spacings so that the cells have 
unitary width; moreover a histogram will be denoted by means of a set of n 2 — 1 
integer {rik} such that 

^n fc = n 2 -l, (23) 
k 

and such that the number of elements per cell nu is the number of spacings be- 
tween k — 1 and k. The elements of a vector A, not necessarily decomposable, 
denote a partition of n 2 — 1; they are also the coordinates of a point in the pos- 
itive hyperoctant of 7?A ~°° which is located in an hypercube whose edges have 
unitary length and whose vertices have integer and nonnegative coordinates. If 
we permute the elements of A, keeping fixed the coordinate system, the new 
vector will be, in general, in another hypercube, similar to the previous one 
and still consistent with the given histogram. It is easy to prove that the total 
number of hypercubes compatible with a given histogram is 

(» 2 - '» ! (24) 



ni!n 2 ! . . . n k \ 

The cube in the positive hyperoctant defined by the relation 

0<Xi<n 2 -l i= 1,2, ...,n 2 - 1 (25) 

contains all the hypercubes compatible with histograms physically meaningful. 
Given a basis of K p and the solutions of inequalities ( |22] ) , one may ask which 
unitary cubes contained in the cube j25| ) intersect K p , i.e. which histograms are 
compatible with such a permutation. Clearly it is enough to consider kernels of 
matrices D p obtained permuting the rows of the matrix (|^) but not its columns: 
permuting the element of a vector A is tantamount as permuting the columns 
of (|5|); therefore if a kernel K p of D p = CT p , where T p is obtained permuting 
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Figure 2: Distribution of spacings of a decomposable non sandwich system. 

only the rows of T, intersects one of the cubes (^4|), so does any K p in which 
the columns of T p have been permuted with an arbitrary permutation. 

Solving inequalities ( 
— in some simple cases we have been able to solve the chain ( |22| ) by means of a 
computer (see fig. ||) — , we would still be faced with the problem of analysing 
the ensemble of histograms so generated. 

5 Conclusions 

The example of the sandwich systems and the general study performed above 
suggest that the weight of decomposable systems whose spacing distributions 
are different from Poisson law is not negligible. The chain of inequalities ( p2| ) 
can be a starting point to give an answer. However our results are sufficient 
evidence that one cannot consider the distribution of spacings of the levels of 
a quantum system as the "signature " of the so called quantum chaos. It is 
worth to note that the we have not considered the semiclassical limit. Even if 
sometimes the behaviour for h — > can give hints on the properties of a system, 
it can also hide its very quantum nature. 
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